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Abstract: We demonstrate that the application of an external magnetic field could lead to an im¬ 
proved background rejection in neutrinoless double-beta (Ovj3j8) decay experiments using a high- 
pressure xenon (HPXe) TPC. HPXe chambers are capable of imaging electron tracks, a feature that 
enhances the separation between signal events (the two electrons emitted in the Ovj3j3 decay of 
^^^Xe) and background events, arising chiefly from single electrons of kinetic energy compatible 
with the end-point of the OvjSjS decay (Qpp). Applying an external magnetic field of sufficiently 
high intensity (in the range of 0.5-1 Tesla for operating pressures in the range of 5-15 atmospheres) 
causes the electrons to produce helical tracks. Assuming the tracks can be properly reconstructed, 
the sign (direction) of curvature can be determined at several points along these tracks, and such in¬ 
formation can be used to separate signal (Ovj8j3) events containing two electrons producing a track 
with two different directions of curvature from background (single-electron) events producing a 
track that should spiral in a single direction. Due to electron multiple scattering, this strategy is not 
perfectly efficient on an event-by-event basis, but a statistical estimator can be constructed which 
can be used to reject background events by one order of magnitude at a moderate cost (~30%) in 
signal efficiency. Combining this estimator with the excellent energy resolution and topological 
signature identification characteristic of the HPXe TPC, it is possible to reach a background rate of 
less than one count per ton-year of exposure. Such a low background rate is an essential feature of 
the next generation of Ovj3j3 experiments, aiming to fully explore the inverse hierarchy of neutrino 
masses. 
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1. Introduction 

Neutrinoless double beta decay (OvjSjS) is a postulated very slow radioactive process in which two 
neutrons inside a nucleus transform into two protons, emitting two electrons. The discovery of this 
process would demonstrate that neutrinos are Majorana particles and that total lepton number is not 
conserved in nature, two findings with far-reaching implications in particle physics and cosmology 
[ 1 , 2 ]. 

After 75 years of experimental effort, no compelling evidence for the existence of Ovj 8 j 8 decay 
has been obtained. The current generation of experiments, with fiducial masses in the range of 
100 kg of isotope and a total background count in the region of interest (ROI) around Qpp of 
a few tens of counts per year, will barely explore the so-called degenerate hierarchy of neutrino 
masses (mi ~ m 2 ~ m 3 ). Exploring the full inverse hierarchy of neutrino masses would require a 
sensitivity to the neutrino Majorana mass of 20 meV, which in turn implies building detectors with 
fiducial masses in the range of one ton of isotope and a total background count in the ROI of, at 
most, a few counts per ton-year of operation. This is a tremendous experimental challenge which 
requires increasing the target mass of the current generation of experiments by typically one order 
of magnitude while at the same time decreasing the backgrounds by at least one order of magnitude 
(for a recent discussion and extensive bibliography on the subject, see [3]). 

One of the technologies currently being developed is that of high pressure xenon (HPXe) 
Time Projection Chambers (TPCs). In particular, the NEXT collaboration [4] is building a HPXe 
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TPC capable of holding 100 kg of xenon enriched at 90% in the j8j8 decaying isotope ^^^Xe. 
NEXT operates at 15 bar and uses electroluminescent (EL) amplification of the ionization signal 
to optimize energy resolution. The detection of EL light provides an energy measurement using 
60 photomultipliers (PMTs) located behind the cathode (the energy plane) as well as tracking via a 
dense array of about 8,000 silicon photomultipliers (SiPMs) located behind the anode (the tracking 
plane). 

The NEXT experiment has completed the R&D phase which was carried out with the large- 
scale prototypes NEXT-DEMO and NEXT-DBDM. The prototypes have measured an energy res¬ 
olution which extrapolates to 0.5-0.7 % EWHM at Qpp- NEXT-DEMO has also shown the robust¬ 
ness of the topological signal [5, 6, 7, 8]. Currently, the NEXT collaboration is commissioning the 
first phase of the experiment, called NEW, at the Canfranc underground laboratory (ESC). NEW is 
a scaled-down version of NEXT-100, deploying 10 kg of isotope mass and 20% of the sensors of 
NEXT-100 (12 PMTs in the energy plane and 1,800 SiPMs in the tracking plane). 

The NEXT background model, which will be fully validated with NEW, predicts a background 
rate of 5 x 10^"^ ctskeV^' kg^^ yr^' in the ROI [9]. The energy resolution for NEXT-100 is 
assumed to be 0.7% EWHM at Qpp (~ 17 keV). NEXT-100 is scheduled to begin acquiring data in 
2017. The experiment expects less than one count of background per 100 kg and year of exposure, 
and thus its sensitivity to is not dominated by background subtraction and increases rapidly 
with exposure. The expected sensitivity to the OvjSjS half-life is 7°,^ > 7 x 10^^ yr for a exposure 
of 300 kg-yr. This translates into a nipp sensitivity range of [90 — 180] meV, depending on the 
nuclear matrix element. An important advantage of the HPXe technology is the fact that it can be 
extrapolated to ton-scale target masses. 

A central feature of a HPXe TPC is the capability of imaging electron tracks, a feature that 
can be used to separate signal events (the two electrons emitted in a 0vj3 j8 decay) from background 
events (mainly due to single electrons with kinetic energy comparable to the end-point of the 0vj3 j3 
decay, Qpp). In this paper, we show that adding an external magnetic field of sufficienlly high 
infensify (in fhe range of 0.5-1 Tesla for operafing pressures in fhe range of 5-15 afmospheres) 
allows for improved separafion of signal and background evenfs provided fhaf fwo key condifions 
are mef. 

1. The spafial resolufion needs fo be sufficienlly good (in fhe range of 2-3 mm). 

2. The elecfron hack needs fo be reconslrucfed wifh sufficienf precision. 

If fhe above fwo condifions are mef, Ihen a slalislical estimator based on fhe average sign of fhe 
curvalure of fhe elecfron hack in fhe magnetic field can be consfrucfed. This estimator can be used 
to rejecf background evenfs by one order of magnilude while mainfaining a 60-80% signal effi¬ 
ciency (depending on fhe pressure). Combining Ibis eslimalor wifh fhe excellenl energy resolution 
and fopological signafure characteristic of fhe HPXe-EL TPCs, we estimate fhaf a fulure Ion-scale 
detector will be able fo reach less fhan one counf background even! per Ion year of exposure, and 
fhus meef fhe requiremenfs needed fo explore fhe inverse hierarchy. 

This paper is organized as follows. Secfion 0 describes fhe reconshucfion of fhe hacks of high- 
energy elecfrons produced in fhe dense gas of a high pressure TPC. The simulafion and analysis 
of such hacks in fhe presence of an exfernal magnelic field is described in secfion ^ In section 
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the analysis procedure is discussed, and the Monte Carlo data sample to which it is applied is 
described in section ^ Finally, the main results are shown in sections ^ and 0, and their implications 
are discussed in section 0. 

2. Imaging tracks in a HPXe TPC 



Figure 1. Principle of operation of an asymmetric HPXe TPC with EL readout. 

Figure |I] shows the principle of operation of an asymmetric HPXe TPC using proportional 
electroluminescent (EL) amplification of the ionization signal (as is the case for NEXT-100). The 
detection process involves the use of the prompt scintillation light (^i) from the gas as the start-of- 
event time, and the drift of the ionization charge to the anode by means of an electric field (~ 0.3 
kV/cm af 15 bar) where secondary EL scinfillafion (^ 2 ) is produced in a narrow region defined 
by fwo highly fransparenf meshes (called fhe EL gap or LEG). High volfages are applied fo fhe 
fwo meshes fo esfablish an elecfric field of ~ 20 kV/cm af 15 bar in fhis region. The defecfion of 
EL lighf provides an energy measuremenf using PMTs in fhe case of NEXT-100 locafed behind 
fhe cafhode (fhe energy plane). The reconsfrucfion of fhe frack fopology is carried ouf wifh a 
dense dense array of SiPMs locafed behind fhe anode (fhe tracking plane). The x-y coordinafes are 
found using fhe informafion provided by fhe fracking plane, while z is defermined by fhe driff time 
befween fhe defecfion of and S 2 . Eor each reconstrucfed space poinf, fhe defecfor also measures 
fhe energy deposifed. Thus, fhe frack is imaged as a collecfion of hifs. Each hif is defined by a 3D 
space coordinafe and by an associafed energy deposifion, as (x,y,z,E). 

Double befa decay evenfs leave a disfincfive topological signafure in a HPXe TPC: a confin- 
uous frack wifh larger energy depositions {blobs) af bofh ends due fo fhe Bragg-like peaks in fhe 
AElAx of fhe sfopping elections (figure 0, leff). In confrasf, background elections are produced by 
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Figure 2. Monte Carlo simulation of a signal (Ovj3j3) event (left) and a background event (right) in xenon 
gas at 15 bar. The color corresponds to energy deposition in the gas. The signal consist of two electrons 
emitted from a common vertex, and thus it features large energy depositions (blobs) at both ends of the track. 
Background events are, typically, single-electron tracks (produced by photoelectric or Compton interactions 
of high energy gammas emitted by or ^os-pj isotopes), and thus feature only one blob. 



Energy Blob 1 Energy Blob 1 


Figure 3. Probability distribution of signal (left) and background (right) events in terms of the energies of 
the end-of-track blobs. The blob labelled as ‘1’ corresponds to the more energetic one, whereas ‘blob T 
corresponds to the less energetic of the two. In a signal event, the blobs have, on average, the same energy. 
In a background event, blob 1 has an energy similar to that of a signal event while the energy of blob 2 is 
very small. 


Compton or photoelectric interactions, and are characterized by a single blob in one of the ends (fig¬ 
ure 0, right). Reconstruction of this topology using the tracking plane provides a powerful means 
of background rejection. For each track, the energy in both extremes of the track is measured and 
labelled as Eb\ (the energy of the most energetic energy deposition), and £"^2 (the least energetic 
energy deposition). In a signal event, E},\ ~ Eb 2 , while for background events Ei,\ » E^i. Figure 
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0 shows how this feature can he used to separate signal from background. 

In this study we compare the single and douhle-electron track signatures in the presence of 
an external magnetic field and quantify the additional ability to reject single-electron events that is 
gained by examining the curvature of the tracks induced by the field. We note that the ability to 
conduct a detailed study of the reconstructed track is highly dependent on the spatial resolution with 
which the track can be reconstructed in the detector and on developing a reconstruction algorithm 
with minimum bias to find the track extremes and track orientation. 

Spacial resolution, in turn, depends on: 

1. Tracking plane design: this includes the pitch of the SiPMs in the tracking plane as well 
as the SiPM response. Indeed, the use of SiPMs with very low dark current and high gain 
allows one to determine the location of an event by weighting the position of each sensor 
with the light recorded, thus improving dramatically the “digital” resolution, which goes as 
~pitch/\/T2- The digital point resolution corresponding to a pitch of 10 mm (NEXT-100) is 
\o/Vn 3 mm. Using weighted information (e.g, local barycenter algorithms), the point 
resolution improves to about 1 mm. 

2. The width of the EL region: The non-zero width of the EL gap adds an extra resolution term 
which goes like w/ \/T2 where w is the width of the grid. Eor the NEXT-DEMO prototype, 
w ~ 5, so this effect would add a resolution term of 1.5 mm. Eor NEXT-100, w ~ 3.5 reducing 
the resolution to 1 mm. 

3. Diffusion of the drifting electron cloud: both longitudinal and transverse diffusion are high 
in pure xenon (of the order of 10 mmJ^/m). On the other hand, work in progress within the 
NEXT collaboration^ suggests that adding small amounts of cooling gases such as CH 4 or 
CE 4 to pure xenon (at the level of 0.1 % of CH 4 or 0.01% of CE 4 reduce both transverse and 
longitudinal diffusion to some 2.0 mmly/va. 

Concerning reconstruction algorithms, the collection of energy depositions (“hits”) that are 
provided by the tracking plane must be organized into an ordered track in which the beginning and 
end of the track are correctly identified. This is a difficult task due to the wide variety of track 
geometries that can result from multiple scattering, and it must be noted that the results of this 
study assume it is done correctly for all events. The NEXT collaboration is currently working on a 
reconstruction algorithm capable of properly ordering the hits with minimum bias towards specific 
geomefries. 

To summarize, it appears possible to achieve a resolution of ~2-3 mm in x-y-z, and to recon¬ 
struct tracks with minimum bias as discussed before. In the remainder of this paper we assume 
both conditions. Our Monte Carlo study smears the individual hits by 2-3 mm and assumes that the 
reconstruction algorithm provides the correct ordering of the track. 

3. Electron tracks propagating in high pressure gas inside a magnetic field 

A particle of charge q moving at a velocity v in the presence of a magnetic field B is acted upon 
by a force F = q{y x B). This force will cause an electron (propagating in vacuum) to execute a 

*“A homeopathic cure to large diffusion in pure xenon”, the NEXT collaboration, paper in preparation 
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helical motion in the magnetic field such that if the thumb of the right hand is positioned along 
the direction of the field line B fhe elecfron will rofafe in fhe direcfion in which fhe fingers curl 
when closed around fhe field line. The frequency of rofafion abouf fhe field line is known as fhe 
cyclofron frequency tUcyc = qB/m, where B is fhe magnifude of fhe magnefic field B = |B| and m is 
fhe mass of fhe charged particle. The radius of curvafure r of fhe frajecfory is direcfly proporfional 
fo fhe election tiansverse momenfum pj and fo fhe inverse of fhe magnifude of fhe magnefic field, 
r = pT/qB. 

If fhe direcfion of fhe applied magnefic field is known, fhe curvafure of fhe frack can be cal- 
culafed fo defermine whefher fhe componenf of fhe election velocity along fhe magnefic field is 
parallel or antiparallel fo fhe field. Assuming fhaf fhe magnefic field is direcfed along fhe z-axis, 
we are inferesfed in fhe curvafure of fhe hack in fhe x-y plane as if progresses in z. The curvafure K 
in fhe x-y plane can be calculated for a frack parameferized by fhe coordinafe z as 


{dx/dz) ■ (d^y/dz^) — (dy/dz) ■ {d^x/dz^) 
r 13/2 

{dx/dzY + {dy/dzY 


(3.1) 


Wifh fhis definifion, an election fraveling in fhe direcfion of fhe magnefic field will rofafe 
around fhe field lines wifh positive curvafure, while an election fraveling opposite fhe direction of 
fhe magnetic field will spiral wifh negafive curvafure. 

The curvafure, however, will be of fhe opposife sign if fhe frack orienfafion is nof property 
idenfified in fhe calculafion (i.e., if dz is of fhe wrong sign). Thus, when calculafing fhe curvafure 
of a single-election frack, one would expecf fc > 0 for dz > 0 and fc < 0 for dz < 0 given fhaf dz 
is always in fhe direction of fhe election velocify. However, for a Ovj3j3 “double election” frack 
(defined as fhe frack befween fhe fwo blobs fhaf mark fhe sfarf and end of fhe frajecfory), faking one 
of fhe extiemes fo be fhe beginning of fhe hack and fhe ofher fo be fhe end will lead fo a calculation 
of K assuming fhe wrong frack orienfafion for one of fhe fwo elections, as fhe verfex af which fhe 
reacfion occurred is found somewhere on fhe inferior of fhe frack. Therefore one expecfs fo find 
ff > 0 for dz < 0 and ff < 0 for dz > 0 for a significanf fracfion of fhe hack. This difference in 
fhe behavior of fhe calculafed curvafure of reconsfrucfed fracks can be used fo disfinguish single- 
elections from 0vj3j3 double elections. 

One serious limifafion fo fhis stiafegy is fhe presence of election mulfiple scaffering (MS). 
This is a process by which an election is scaffered repeafedly while fraveling fhrough fhe dense gas, 
resulfing in deviafions of ifs pafh by significanf angles. The process can be described approximafely 
by considering fhe angle of deflection projecfed on fhe x and y planes when scaffered fhrough a 
fhickness of xenon x (assuming fhaf fhe election is fraveling along an axis in fhe z direcfion). Each 
of fhese angles Qx and dy can be modeled by a gaussian disfribufion (Moliere fheory) as 


o'^{ex,y) = \/dz/Lo [1 + 0.038ln(dz/Lo)], (3.2) 

for an election wifh momenfum p in MeV/c before scaffering and befa factor j3 (where c = 1). Lq 
is fhe radiafion lengfh, a property of fhe medium. 

In fhe absence of MS fhe tiansverse momenfum of an elecfron could be direcfly defermined 
by fhe hack curvafure. Consequenfly, a single-election frack—for which fhe momenfum is greafesf 
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at one end of the track and decreases toward the other end—would be clearly distinguishable from 
a double-electron track consisting of two electrons traveling in opposite directions originating at 
some vertex in the middle of the track. Unfortunately, in a HPXe TPC multiple scattering is high, 
resulting in large errors that complicate the direct measurement of the magnitude of the curvature. 
The sign of the curvature, however, is much less affected by MS. Consequently we choose the sign 
of the curvature, rather than its magnitude, as our discriminating variable. 

4. Determination of the Track Curvature 

Given a track of electron ionization defined as a collection of hits containing position and energy 
information, {x,y,z,E), the curvature at each hit in the track is calculated by numerical computation 
of the derivatives dx/dz, dy/dz, d^x/dz?', and d^y/dz}. From these the curvature K is calculated at 
each point. Since we do not have x and y as a function of z but rather x, y, and z as a function of hit 
number n, we can calculate the derivatives x' = dx/dn, y' = dy/dn, z! = dz/dn using the chain rule 
as ^ = x' /z!, and 

d^x x!' —z!'{dx/dz) 

The expressions for dy/dz and d^y/dz^ can be obtained by replacing in the above x —)• y. 
Note that outliers may need to be removed from the resulting arrays of first and second derivatives 
due to points between which the z-coordinate changes very little. To ensure more stable values 
of the derivatives, an outlier removal procedure is applied to all derivatives and second derivatives 
computed which consists of iteratively calculating the mean and variance a of each array, replacing 
any value that lies outside of 5a of the mean value with the average of the two nearest values in the 
array, and continuing this procedure until the calculated variance is no longer less than the value of 
the variance calculated in the previous iteration. 

The curvature calculated using each pair of points is then corrected as follows: if for the two 
points Z 2 < Zi, that is dz < 0, the curvature is multiplied by -1 (see section Note that the outlier 
removal procedure described above is also applied to the calculated curvature array. 

The discriminant variable is obtained by condensing the curvature information for each track 
into a single numerical value indicative of the nature of the track. First a curvature sign array is 
created consisting of values of either +1 or — 1 depending on the sign of each value in the calculated 
curvature array. The numerical value, which we will call the curvature asymmetry factor (CAF), is 
defined as the average of the curvature sign array using elements in the first half of the track minus 
the average of the curvature sign array using elements in the second half of the track: 

1 /iV/2-1 N \ 

^C = TJ^{ £ sgn(Kjt)- £ sgn(K-i) . (4.2) 

\ k=0 k=N/2 } 

Note that if the number of hits in the track is odd, the first term includes the first (N — 1) /2 hits and 
the second term includes the remaining (N—l)/2 + l hits. Though due to multiple scattering and 
imperfect reconstruction resolution, this factor is not a perfect discriminator, statistically it can be 
used to separate single-electron and double-beta events. 
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To further illustrate this idea, we consider a simple toy Monte Carlo in which tracks are gen¬ 
erated in steps of 10 keV energy deposition using the NIST [10] dE/dx curve and equation 
to model multiple scattering. For each energy deposition AF = 10 keV, the electron is prop¬ 
agated along a path, calculated according to the force exerted by the magnetic field, of length 
Ax = {dE/dx)^^dE. The initial direction of travel is assigned at random and is updated af¬ 
ter each step Ax by calculating a random dx and 6y using equation and assuming the direction of 
travel is along the z-axis. A single-electron track is generated in this way assuming an initial energy 
of 2.447 MeV, to simulate background events. For “signal” events, two single-electron tracks are 
generated from a common vertex, each with an initial energy of 2.447/2 = 1.2235 MeV. Notice 
that this is not realistic, since the two electrons produced in Ovj3j8 events share their energy only 
statistically. Event-by-event, each electron can have an energy that varies between 0 and Qpp with 
the other electron taking the rest. This makes it much harder, in practice, to separate Ovj3j8 events 
from single electrons than in the naive case assumed here for illustration purposes. 

Our toy Monte Carlo allows us to study the influence of (fhe gaussian componenf of) multiple 
scattering presenf in fhe simulafed evenfs in a sfraighfforward manner by changing fhe pre-facfor 
in equation from 13.6 MeV fo a lower or higher value. Figure ^ shows several signal and 
background fracks and fheir calculafed curvafure sign simulafed wifh differenf mulfiple scatter¬ 
ing pre-facfors corresponding fo no mulfiple scattering (0 MeV), fhe sfandard amounf of multiple 
scattering (13.6 MeV) and fwice fhe sfandard amounf of mulfiple scattering (27.2 MeV). 

Increasing mulfiple scattering makes fhe background rejecfion procedure more difficulf, as 
fhe random motion of MS offen works againsf fhe curvafure infroduced by fhe magnefic field and 
quanlified in equation However, as shown in figure 0 fhe CAF separafes single-election and 
double-election fracks perfecfly, even when fhe amounf of MS is doubled. Nafurally, fhis resulf is 
too opfimisfic, as already mentioned, for several reasons: fhe foy MC only includes fhe gaussian 
componenf of MS, neglecfing imporfanf non-gaussian effecfs (included in fhe Geanf-4 simulafion), 
fhe poinf resolufion is assumed perfecf, and fhe "double-elections" share exacfly half of fhe energy. 
However, fhe good separation shown in figure ^ demonstrates that the proposed procedure is sound. 

5. Data sample 

The primary simulation-based study consisted of analysis of Monte Carlo datasets generated in a 
large virtual box of high pressure xenon gas using GEANT4 [11]. Such an arrangement emulates 
events contained in the fiducial volume of NEXT-100. Eor various configuralions of gas pressure 
and magnefic field, 10^ signal (OVj8j3) evenfs and lO'* background evenfs (single elections of kinefic 
energy equal to Qpp) were generafed. 

Each simulated frack was recorded as a series of hifs consisting of a location (x,y,z) and a 
deposifed energy E. Space coordinafes were grouped in a single continuous frack consfrucfed by 
adding fogefher all fhe hifs leff by fhe election (in fhe case of single election evenfs) or fhe fwo 
elections (in fhe case of 0Vj8j3 evenfs). This procedure emulates fhe pattern recognition in a HPXe 
TPC, which reconsttucfs as a single frack bofh single and double elections (see figure @), and 
assumes fhaf fhe hifs are reconsfructed in fhe proper order. Eurfhermore, fhe hifs were smeared fo 
simulafe poinf resolufion. Before proceeding wifh furfher analysis, fhe summed energy of all hifs 



Figure 4. Calculated curvature sign at each point along the track for a single-electron event (left) and a 
(simplified, see text) double-beta event (right). The tracks were generated with no multiple scattering (top), 
the standard amount of multiple scattering (middle) and twice the standard amount of multiple scattering 
(bottom). The red + markers indicate positive curvature, while the blue dots indicate negative curvature. 


recorded in the track was required to be at least 2.4 MeV. This cut corresponds to demanding the 
event energy to be within the ROI. 

The X, y, and z-coordinates of the track were then separated into individual arrays, and a low- 
pass filter was applied to each array. This served to smooth the track and facilitates the calculation 
of the derivatives. Figure ^ shows an example of a track before and after the filtering procedure, 
the details of which are discussed in appendix ^ 
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Curvature asymmetry factor Curvature asymmetry factor 4>c 


Figure 5. CAF calculated for 10^ single-electron and double-electron tracks generated with the standard 
amount of multiple scattering (left) and twice the standard amount of multiple scattering (right). Note that 
though the background rejection capability appears almost perfect based on the results of this simple Monte 
Carlo, other factors such as the reconstruction resolution and non-equal energy sharing between the two 
electrons in Ov)3j3 decay degrade the discrimination power and are not taken into account here. 



Figure 6. An example of a single-electron track before (left) and after (right) the application of a lowpass 
filter. The filter applied to this track is shown in appendix ^ figure |T^ 


6. Results 

The analysis described in section ^ was performed for single-electron and Ovj3j3 tracks for all 
combinations of P = 5, 10, and 15 atm and B = 0.1, 0.3, 0.5, 0.7, and 1.0 T. In all of these 
configurations, the hits were smeared in {x,y,z) with Os = 2 mm. They were also “sparsed,” that 
is, groups of every Ns hits were averaged into a single hit, with Ns = 2. Additional analyses in the 
P = 10 atm, B = 0.5 T configuration were performed, one with = 1 mm and Ns = 1, and another 
with Os = 2 mm and Ns = 3. For each configuration, 10^ events of each track type (single-electron 
and 0vj3j8) were generated. For illustrative purposes, figure 0 shows the calculated curvature signs 
for a single-electron and double-beta event, after filtering with a lowpass filter, at P = 0.5 T and 
P = 10 atm. 

For each track we calculate the curvature asymmetry factor shown in equation Next, a 
cut was placed on the CAF, defining the events considered to be candidate OvjSjS tracks. The cut 
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Figure 7. Calculated curvature sign at each point along the track for a single-electron event (left) and 
a double-beta event (right) generated using the GEANT-4 based Monte Carlo. The red -f markers indicate 
positive curvature, while the blue dots indicate negative curvature. Note that a lowpass filter has been applied 
to the tracks shown as described in appendix ^ 

was varied, and in each case the fraction of single-electron (background) tracks rejected by the cut 
was determined along with the fraction of 0Vj8j3 (signal) tracks accepted.^ An example of such an 
analysis is summarized in figure H. 




Figure 8 . Curvature asymmetry factor for single-electron and OvjSjS events with (7j = 2 mm (left) and 
resulting signal efficiency vs. background rejection curve produced by varying a cut on 0c space, shown 
also for (7j = 1 mm and = 3 mm (right). Here s is the fraction of the total number of Ovj3j3 events 
considered that was identified as a OvjSjS event, and b is the fraction of the total number of single-electron 
background events considered that was identified as a Ov)3)3 event. These results were obtained for the 
/* = 10 atm and B ~ 0.5 T configuration. 

To illustrate the performance of the method at the different configurations of gas pressure and 
magnetic field, we examine the signal efficiency s (fraction of OVj8j3 events correctly identified) 

^Of the 10^ events generated for each configuration and track type, only events that passed a set of initial topological 
cuts were considered in the curvature analysis. These cuts required a minimum total deposited energy of 2.4 MeV and 
the deposition of this energy in a single continuous track, so events with some energy deposited significantly away (> 2 
cm) from the main track were not analyzed further. The final fraction of events therefore refers to a fraction of this subset 
of events analyzed, not a fraction of the total 10^ events generated. 
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Figure 9. Signal efficiency corresponding to approximately 80% background rejection (above) and 90% 
background rejection (below) vs. magnitude of applied external magnetic field for different pressures. 


obtained given a cut that provides approximately 80% and 90% background rejection. The results 
are shown in figure ^ 

Note that the results depend significantly of the operating pressure. The best performance is 
obtained, as expected, at the lowest considered pressure (5 bar), and it reaches a plateau (80% effi¬ 
ciency for a background rejection of 80%, and 70% efficiency for a background rejection of 90%) 
for magnetic fields larger than 0.5 Tesla. The performance at 10 bar appears also quite accept¬ 
able. It reaches a plateau (75% efficiency for a background rejection of 80%, and 65% efficiency 
for a background rejection of 90%) for magnetic fields larger than 0.7 Tesla. The performance 
at 15 bar is poorer and does not saturate with the magnetic field, suggesting that some additional 
improvement can be obtained with larger values of B. 
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7. An improved analysis 


While the curvature asymmetry factor (^c from equation provides a reasonable separation be¬ 
tween signal and background events, it still represents a rather minimal approach. For example, in 
calculating (j)c, the difference was taken in the average of the signs over the two halves of the track, 
but the vertex in a OVj8j3 event is not likely to be found in the center. A more detailed analysis 
taking into account more features of the two different classes of events should give improved per¬ 
formance. As an example of such an analysis, we consider constructing a track “profile” for each 
type of event containing the value of the curvature sign at each point along the track averaged over 
all events. Since not all tracks will be of exactly the same length, we normalize the track length to 
1 , and thus the profile confains an average sign value for each fractional disfance along fhe frack. 
Nofe fhaf fraclional disfance is defermined by hif number k divided by fhe fofal number of hits N 
rather than actual distance traveled. Figure [T^ shows profiles made for simulated single-elecfron 
and Ovj3j8 evenfs. The profiles show a clear separafion befween signal and background evenfs 
earlier in fhe frack where fhe sign of fhe curvafure is mosf likely fo differ. 



Figure 10. Average curvature sign profiles for single-electron and Ov)3)3 events generated at B = 0.5 T and 
P = 10 atm, with Oi = 2 mm and W = 2. The profile shows the average sign vs. fractional distance along the 
track {k/N, where k is the hit number and A is the total number of hits). The profiles were created by dividing 
the normalized track into 15 bins, and for each calculated curvature value for each event, adding a value of 
-1-1 or —1 to the bin, then dividing each bin by the total number of values placed in it. The values in the first 
and second bin were extrapolated linearly to k/N = 0, and those in the final two bins were extrapolated to 
k/N = 1. The entire set of values was then interpolated using cubic polynomials. 


Once the profiles have been generaled, Ihey can be used lo define a new variable lhal will 
provide separation befween fhe Iwo classes of evenfs. For each evenl, fhe sign of fhe curvafure 
sgn(K';) af each (normalized) disfance along fhe frack is compared fo fhaf of fhe profile for signal 
(OVj8j3) evenfs, ps, and background (single-elecfron) evenfs, ps, in creating fwo ;^^-like sums. 


2 ^ y -Psjk/N)]'^ 

to 


(7.1) 
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and a similar sum Xb which the background profile ps is used in place of ps- The values ps{x) 
and pb{x) are the mean curvature signs at fractional track location = k/N, and the values a| 
and aj are the variances of the mean curvature sign. The variances were found to be near 1 at 
all k/N in both profiles, reflecting the large amount of fluctuation in curvature sign at any given 
point on the track due to multiple scattering. Despite these large variances, the large number of 
events used in determining the average values made for low errors on the mean and thus smooth 
profiles. The discriminating variable is then = Xs/Xs- background events, this ratio should 
be smaller, as the curvature signs throughout the track should be more similar to the background 
event profile (smaller Xb) more different than the signal event profile (larger Xs)- Based on 
similar arguments, this ratio should be larger for signal events. 

Figure compares the signal efficiency obtained for a given background rejection using the 
profile-based method (discriminant variable described in this section and the method based on 
the asymmetry factor described in section The profile-based method performs better up to 
a background rejection of about 95%. This clearly shows the potential for further optimization. 
In our “standard test case” (10 atm, 2 mm resolution, 0.5 Tesla) a rejection factor of one order of 
magnitude is achieved for a signal efficiency of about 75% using the profile method. 



Figure 11. Signal acceptance efficiency vs. background rejection obtained at B = 0.5 T and P = 10 atm, 
with (Jj = 2 mm and = 2. The curve obtained using the curvature asymmetry factor is shown along with 
that obtained using the profile-based method. 


8. Outlook: the case of a ton-scale HPXe TPC with a magnetic field 

The application of an external magnetic field in a high-pressure detector capable of accurate particle 
track reconstruction with resolution of approximately 2 mm in (x,y,z) can provide an additional 
background rejection factor of 90% with good signal efficiency (70-80%) at pressures of about 
10 atm. The background rejection improves with decreasing pressure due to the more extended 
tracks produced at lower gas densities. A sufficient track reconstruction algorithm has not yet been 
realized to take full advantage of these findings, and the NEXT collaboration is currently working 
towards a potentially suitable algorithm. 
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It follows that the background rejection of a future ton-scale HPXe using a magnetic field 
could be up to a factor of 10 better than the current rejection factor expected by the NEXT-100 
experiment (while keeping an efficiency in fhe range of 25%). The resulfing background rafe 
would be around 5 • 10^^ cfs keV^^ kg^' yr^^ or abouf a half counf per ton and year in fhe ROI, 
assuming a resolufion of 0.5% FWHM af Qpp- This is enough, for mosf NME sefs, to fully cover 
fhe inverse hierarchy in a few years. 

A. A Lowpass FIR Filter 

One way fo make fhe calculation of fhe curvafure of a frack less suscepfible fo noise infroduced 
fhrough mulfiple scaffering and non-ideal reconsfrucfion resolution is fo apply a lowpass filler fo 
fhe lisfs of values for each coordinale x, y, and z. We now look al fhe arrays of x, y, and z coordinales 
as digilal signals in fhe lime domain, e.g. x[n], where n is fhe index of fhe corresponding hil. Each 
array can be represenled in fhe frequency domain X[k] using fhe discrete Eourier Iransform 

N-l 

X[k] = £ (A.1) 

n—0 

where N is fhe folal number of samples and k is fhe discrele frequency of each complex sinusoid 
^-iinkn/N cycles per A samples. This discrele frequency can be Iranslaled to an analog frequency 
fk (for example in unifs of lime^^) by knowing fhe frequency al which Ihe digilal signal was 
sampled, or Ihe sampling frequency fs in samples per unil time, as 

fk = kfs. (A.2) 

The goal is to apply a digilal lowpass filler to Ihe coordinate arrays lhal serves to smoolh Ihe 
Irack by eliminating high-frequency noise yel relains Ihe curvalure of Ihe Irack due to Ihe magnetic 
field. The ideal lowpass filter will serve to eliminate sinusoidal componenls X[k] in Ihe digilal 
signal for k greater lhan some cutoff frequency kc and allow olhers wilh k less lhan kc- In practice, 
Ihe filler will have a non-ideal stopband in which sinusoidal componenls wilh frequencies less 
lhan kc will be increasingly preserved wilh decreasing k and Ihose wilh frequencies greater lhan kc 
will be increasingly attenuated wilh increasing k. The filter musl be designed to ensure we donT 
eliminate Ihe sinusoidal motion infroduced by Ihe magnetic field, and we use a filler wilh a wide 
stopband (Ihis reduces ils complexify) and place kc near Ihe discrete frequency corresponding to 
Ihe cyclolron frequency 


tocyc = qB/nie = 1.76B x 10^^ rad/s, (A.3) 

where q k, —1.60 x 10^^^ is Ihe eleclron charge in Coulombs, B is Ihe magnetic field slrenglh in 
Tesla, and ps 9.11 x 10^^^ is Ihe eleclron mass in kg. The filler will be a finite impulse response 
(EIR) filler and will lake Ihe form of a lisl of coeflicienls bm- It is applied to the signal x[n] as 


N, 

Xf[n] = ^ bmx[n — m\ (A.4) 

m=0 
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where the number of filter coefficients Nt is the order of the filter. A delay will exist in the final 
filfered signal equal to Nt/2 samples, and fhe final Nt/2 samples will nof be useful in fhe final anal¬ 
ysis and are removed. For a frack sampled in x, y, and z for a fofal of N samples, fo defermine whaf 
discrefe frequency kcyc to which this analog frequency (C0cyc/27r in cycles per second) corresponds, 
one must know the frequency with which the track has been sampled. An average track production 
time T can be calculated from tabulated values of dE/dx (which, since this refers to kinetic energy 
loss, we will call dK/dx) as 


T 


1 

c Jo 


K + me 


{dK/dx) 


dK 1.25 ns, 


(A.5) 


where here we have used Kf = = 2.447 MeV and tabulated dE/dx from NIST [10] in xenon. 

The sampling frequency is then fs = N/T, and so the motion due to the magnetic field should 
manifesl ilself in fhe arrays of sampled x and y coordinafes of fhe frack as a sinusoidal componenf 
of discrefe frequency 


^cyc — {qB/m)-{T/N). (A.6) 

One should ensure fhaf fhe frack is sampled af a rafe higher fhan fwo limes fhe cyclolron frequency, 
lhal is N/T > (Ocyc/2n, or fhe helical mofion will nof be properly reconsfrucfed. An example filler 
is shown in figure |T^ 



Figure 12. FIR filter gain vs. discrete frequency k (solid black curve) designed for a track of length A = 198 
samples. The coefficients X\k\ of the fast fourier transform (FFT) of the x-coordinate array of the track are 
shown vs. discrete frequency in red. The cutoff frequency used in designing the filter is shown as a blue 
vertical line. The filter serves to eliminate the high-frequency noise components of the track that make the 
calculation of derivatives more difficult. 

For the track shown in figure an FIR filter was designed with cutoff frequency equal to 
kc = 1.2kcyc and stopband width equal to 0.2. This was sufficient to eliminate high-frequency 
noise, leaving a smooth track for the calculation of derivatives, yet the filter order was low enough 
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so as to not lose a significant number of the track samples. For this track, a number of samples 


with B = 0.5 T, we have = 0.089 cycles/sample as shown in figure 

As fhe magnifude of B increases, fhe frequency of helical motion oflen becomes greaf enough 
fhaf placing fhe filler culoff frequency kc near fhe cyclolron frequency does nof permil sufficienl til¬ 
lering of fhe high-frequency componenls of fhe Irack, and fhe resulling Irack is nof smoolh enough 
for a clean calculalion of fhe curvalure. Therefore if may be an advanlage fo place fhe culoff fre¬ 
quency of fhe filler below k^yc despife increased allenualion al fhe discrefe frequency correspond¬ 
ing lo fhe helical motion induced by fhe magnelic field.^ Figure |I^ shows fhe signal efficiency vs. 
magnelic field for a background rejeclion of 90% al 10 aim in pure xenon for which fhe filler was 
designed for each Irack individually wilh culoff frequency equal lo 1 .2kcyc, and for which Ihe culoff 
frequency of all fillers was chosen lo be a fixed value. 


A = 198 and using T = 1.25 ns from equation 


aiiva u^c 





Figure 13. Signal efficiency vs. magnetic field for 90% background rejection shown with two different filter 
calculations to demonstrate the effect of the lowpass filter design on the background rejection performance. 
The use of a filter with cutoff frequency calculated on a track-by-track basis (“variable filter”) does not yield 
a remarkable improvement in performance at low values of B and allows too much noise to remain at high 
values of B. The use of a filter with a fixed cutoff frequency (“fixed filter,” same as Xe 10 atm curve in figure 
^ left) in the range of those calculated for B = 0.5 T for all tracks eliminates the noise but does not eliminate 
information about the curvature to the extent that the performance is degraded. 


Because the fixed-cutoff analysis performs significantly better at higher fields, if was chosen 
as fhe primary analysis in fhis sfudy. For fhe resulfs shown in section fhe lowpass filler applied 
lo all Iracks was designed wilh a fixed culoff frequency of kc/ln = 0.102 (wilh fhe exceplion of 
Ihose shown in figure 0, lo which fhe filler shown in figure [T^ specifically was applied). 


^In any case, due to the constant diversion from the ideal path of the electron introduced by multiple scattering, the 
curvature of the track in the magnetic field will not appear at a single distinct value in the frequency domain. 
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